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A universality ansatz for the contact values of a multicomponent mixture of additive hard 
spheres is used to propose new formulae for the case of disparate-size binary mixtures. A 
comparison with simulation data and with a recent proposal by Alawneh and Henderson for 
binary mixtures shows reasonably good agreement with the predictions for the contact values 
of the large-large radial distribution functions. A discussion on the usefulness and limitations 
of the new proposals is also presented. 
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1. Introduction 

A well known result in statistical mechanics is the close relationship between the 
structural and thermodynamic properties of fluids. In the particular case of a mul- 
ticomponent mixture of (additive) hard spheres (HS), the compressibility factor 
Z = p/pksT (where p, p and T are the pressure, density and temperature of the 
system and /c^ is the Boltzmann constant) is given by 



Z =1 + 4:r]^-^XiXjgij{aij), (1) 

where is the number of components in the mixture, Xi is the mole fraction 
of species i, aij = ^{ai + aj), ai being the diameter of a sphere of species i, 

T] = (7r/6)p(o"^) is the packing fraction (with (a^) = YliLi^i^?) diji'^ij) ^-^e 
the contact values of the radial distribution functions (RDFs) gij{r). Hence, in 
this system, knowledge of the contact values of the RDFs suffices to completely 
determine the equation of state (EOS). Exact results for gij{aij) are unfortunately 
only known to first order in the packing fraction, where one simply has 
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where 



(3) 



Therefore, in general one must rely on either approximations or values obtained 
from computer simulations. Amongst the most important approximate analytical 
expressions for the contact values are the ones that follow from the solution of the 
Percus-Yevick (PY) equation of additive HS mixtures by Lebowitz [T], namely 
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and the results obtained from the Scaled Particle Theory (SPT) [2, SjH^IHl, 
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However, neither the PY nor the SPT lead to particularly accurate values and so 
Boublik and, independently, Grundke and Henderson [7] and Lee and Levesque 
@] proposed an interpolation between the PY and the SPT contact values, that we 
will refer to as the BGHLL values: 



BGHLL 
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These lead, through Eq. ([T]), to the popular and widely used Boubhk-Mansoori- 
Carnahan-Starling-Leland (BMCSL) EOS ^ for HS mixtures. For binary mix- 
tures the BGHLL values are rather accurate except for asymmetric mixtures in 
which large hard spheres are present in extreme dilution in a small hard-sphere sol- 
vent. Since such mixtures represent good models for colloidal suspensions and there 
is growing interest in these kinds of systems, there have been man y ef forts to correct 
the deficiencies of the BGHLL values [lj,[iiyi3,|Ti,[li,[il,H 

In this paper the recent proposal of Ref. [2^ is of particular concern, so we also 
quote it explicitly here. It relies heavily on the proposal of Viduna and Smith (VS) 
18l | which reads 
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Alawneh and Henderson (AH) [22[ restrict themselves to binary mixtures (A^ = 2, 
o'2 > ci) and propose to keep VS's expression for (711 (cJi) but amend 5'i2(<7i2) and 
522(0-2) as 



9l2 (0^12) = 9l2 ('^12 
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Here R = (J^jcyi [23*]. With these, they get good agreement for the large-small 
contact values and reasonably good agreement for the large-large contact values as 
compared to Molecular Dynamics results for mixtures with size ratios going from 
= 1 to i? = 10. 

lie 



27l | three of us have advocated 



On a different vein, in recent papers [2J, |2£ 
a "universality" approach as a guide to propose simple and accurate approximate 
expressions for the contact values gijipij) of additive HS mixtures. In this approach, 
if the packing fraction r] is fixed, the expressions for the contact values of the RDF 
for all pairs ij of like and unlike species depend on the diameters of both species and 
on the composition of the mixture only through the single dimensionless parameter 



Zij [cf. Eq. ([3])], irrespective of the number of components in the mixture. Note that 
this feature is shared by the PY, SPT and BGHLL results but not by the VS or the 
AH approximations. The aim of this paper is two-fold. On the one hand, we will 
assess whether the new simulation data of Ref. [22] comply with the universality 
ansatz. On the other hand, we will introduce a refinement of some of the simple 



proposals of Refs. [2^, |25|, |26|, |27|] that retain the universality feature but are able 



to improve the performance in the case of highly asymmetric mixtures. 

The paper is organized as follows. In Sect. [2] and in order to make the paper self- 
contained, we provide a rather brief account of the simple proposals based on the 
universality ansatz and a simple extension to deal with highly asymmetric systems. 
This is followed in Sect. [3] by a comparison between the AH contact values, those 
of the new proposals and the simulation data. The paper is closed in Sect. [4] with 
further discussion and some concluding remarks. 



2. Proposals for contact values including the "universality" assumption 

We start by recal ling the main aspects of the "universality" approach that has been 



used recently [2J, |25|, l26|, |27[ as a guide to propose simple and accurate approximate 
expressions for the contact values gijipij) of additive HS mixtures. According to 
this approach, at a given total packing fraction 77, gij{crij) depends on the size and 
concentration distributions only through the single parameter Zij: 

gijiaij) = G{'n,Zij), (10) 

where the function G{ri,z) is a common function for all the pairs regardless 
of the composition and number of components of the mixture. Note that the exact 
contact values to first order in density [cf. Eq. ^] are consistent with the ansatz 
^ with G{rj,z) ^ l + r/(l + 

In order to propose specific forms for G{r], z) use is made of three basic consis- 
tency conditions. First, in the limit in which one of the species, say i, is made of 
point particles (i.e., ctj ^ 0), gaiai) must take the ideal gas value, except that one 
has to take into account that the available volume fraction is 1 — r/. Thus, 



lim quiai) = . (11) 

^,^0^^ ' l-r] ^ ' 

Next, if all the species have the same size, {(Jk} cr, the system becomes equivalent 
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to a single component system, so that 



lim gijicTij) = Qpir]), 



(12) 



where gp{ri) is the contact value of the RDF of the single component fluid at the 
same packing fraction r] as that of the mixture. 

The third condition is the equality between the pressure in the bulk and near a 
hard wall 0,113]: 



id 3 



(13) 



where g^j is the wall-particle correlation function at contact. Since a hard wall can 
be seen as a sphere of infinite diameter |28|], the contact value g^j can be obtained 
from gijiaij) as 



gwj = Jiin gij{(7ij] 



(14) 



In the special case of a pure fluid plus a hard wall, Eq. (jT3|l becomes 



jjin 512(0-12) = 1 +^r]gp{r]) = Zp{r]), 



(15) 



where Zp[r]) is the compressibility factor of the pure fluid. 

Assuming the validity of the ansatz ([TO]) , conditions (fTTj) . (fT2]) . and (fT5 
become 
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respectively. Conditions p^. p7|) and ([T9|) are fulfllled by the quadratic proposal 



1 — r] 
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However, Ge2 does not satisfy Eq. (jlSp for an arbitrary number of components 
A'' > 3. This is accomplished by a cubic form H, 27]: 



Ge3(r?, z) = Ge2{v, z) + {1 - T^) [gf^[r,) - gpir,)] z{z - l)(z - 2), (21) 
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Figure 1. (Colour online) Plot of Z22{a;2) as a function of l/R for, from top to bottom, X2 = 0.001, 0.002, 
0.005, 0.01, 0.02 and 0.05. 



where 



9^ 



1 -r?/2 + 7?V4 
(1-^)3 



(22) 



is the SPT contact value of the pure fluid. The proposals Ge2 and Ges are very 
flexible in the sense that they can accommodate any reasonable choice for gp{ri). 
In particular if gp{'q) = g^^{r]) then Ge2{f],z) = Gez{f],z) and one recovers the 
SPT contact values given in Eq. ([5]). 

Taking either the gp{ri) values that follow from the Carnahan-Starling or the 
Carnahan-Starling-Kolafa (CSK) EOS as input, both approximations Ge2 and 
Ge3 yield in general very good predictions for OjAojj) in binary and polydisperse 
mixtures under a wide range of conditions 24 , 23, 2^, [13] , which typically implies 



Zij < 2. However, in the case of very disparate mixtures, the parameter zij can 
take larger values. For instance, in a binary mixture with i? > 1 one has 



Z22{X2,R) = R 



I- X2+ X2R^ 
1 - X2+ X2R^ ' 



(23) 



The behaviour of 222 as a function of l/R is shown in Fig. [T]for different values 
of X2- It follows that, at fixed X2, there exists a value R = Rmsi.x{x2) at which Z22 
reaches a maximum z^^^{x2)- The peak z^^{x2) becomes narrower and higher 
as X2 decreases, its location moving to the left (larger values of R). An explicit 
expression for Rma.x{x2) can be derived by differentiation of 2:22 with respect to i?, 
which yields a cubic equation for R whose real solution is 



Rmiix{x2) 



1 

2 L 



1 + (2x2 ^ + 2x^Vl -a;2 - 1) + {2x2^ + 2x2 ^\/T 



X2-1 



-1/31 



(24)" 

A plot of Rma.x{x2) and the associated values of Z22'^{x2) is shown in Fig. [2j In 
the limit X2 ^ one gets i?max {2x2)~^^^, so that limaj^^o 3;2-Rmax = 1/2- For 
X2 > one has X2i?max > 1/2- This means that, when Z22 reaches its maximum 
value, the partial volume fraction occupied by the big spheres is always larger than 
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Figure 2. (Colour online) Plot of -Rmax(a;2) (solid line) and the corresponding z^^{x2) (dashed line). 
Also shown is the behaviour of the quantity X2i?max (dotted line, right vertical scale) as a function of X2. 

half the one occupied by the small ones. Figure [2] shows that z^^^{x2) can reach 
rather high values for sufficiently small X2- In particular, z^^{x2) > 2 if X2 < 0.036, 
irrespective of the value of R. Since G{t], z) is a monotonically increasing function 
of z, the larger Z22, the higher 522 (o"2)- Thus, Eq. gives the value of R at 
which (?22(c2) has a maximum, according to the ansatz ([TO]) . Note that i?™^^(x2) 
is independent of rj. 

The polynomial forms (|2U|) and (|21|) are constructed by imposing exact properties 
at z = 0, z = 1 and z = 2. On the other hand, in disparate mixtures where one 
can have values of z higher than 2, Eqs. I\20h and ()2ip are not accurate enough. 
In fact, again in the binary case, there exists evidence indicating that 522 (o'2) 
grows exponentially for large R, as explicitly incorporated in the AH proposal [22], 
Eq. ([9]). This translates, in the context of the ansatz (jlOp . into an exponential 
dependence of G{r], z) on z for large z. Thus, here we propose to amend Eqs. (|20]) 
and (EH) as 



Ge2+ir],z) 



Ge2ir],z), 



0< z<2, 



Zp{r])exp[Ae2ir])iz-2)], 2 < z, 
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Ge3ir],z), 0<z<2, 
Zpirj)exp[Ae3i7]){z-2)], 2 < z, 



(26) 



where Ae2 (f?) and Ags (77) are determined by requiring continuity of the first deriva- 
tive at the matching point z = 2. Their expressions are 



2 - 37]/2 
1 — r] 



2(1 - 3r/)(7p(r?) 
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Note that the amended proposals (j25p and (I26p also share with Eqs. (I20p and 
(|2ip the flexibility to accommodate any reasonable choice for gp. In what follows 
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dash-dot line: Gesfr;, z); dotted line: Ge2+('7, -z); solid line: Ge3+{»?, z); solid circles: simulation data of 
(;il((7i) from Ref. |22j : solid triangles: simulation data of (;i2(o'i2) from Ref. [22;] ; solid squares: simulation 
data of 922 ('T2) from Ref. [22l |: open circles: simulation data of g^i and 3,1,2 from Ref. [29l |. The inset shows 
the region < Zij < 2. Note that Ge2+('?i and Ge3{r], z) are practically indistinguishable at this density. 




we will use the one corresponding to the CSK EOS, namely 



CSK, , _ l-r?/2 + r?2(l- 27?)/12 



,f ^(,) = - -^/--^v^^ (29) 



3. Comparison with simulation data 



Once the new universal proposals have been introduced, the question now arises 
as to whether the available simulation data are consistent with the universality 
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Figure 5. (Colour online) Same as in Fig. |3] but for rj = 0.4. 




Figure 6. (Colour online) Same as in Fig. [S] but for r] = 0.5. At this density there are no simulation data 

for Qyjj. 



ansatz. Were this to be the case, for a given packing fraction rj plots of gij vs Zij 
of different mixtures should lie on a common curve and this will be tested next. 

The approximations Ge2(^, z), Gesif], Ge2+{Vi ^) G^z+iVi given in Eqs. 
(f20]) (f2T]) . ([25]) and ([26ll . respectively, complemented with the CSK expression for 



29(] in 



dpiil) given in Eq. ([29]), are compared with simulation data of 
Figs. [SHSl We can observe that Ge2 and GeS perform very well for < z < 2, 
except for the higher values of 922 (f 2) in that region (see the insets of Figs. [3H6]). 
Since, by construction, both Ge2 and Ge3 give the correct value at z = 2 [cf. Eq. 
Iligp ]. the deviations of the simulation values of 522 (o'2) from the theoretical curves 
in the region < z < 2 can be interpreted as a manifestation of the approximate 
character of the universality ansatz (llOp . This seems to be confirmed by the scatter 
of the simulation data of 522 (o'2) for z > 2. On the other hand, the overall behaviour 
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Figure 7. (Colour online) Plot of 922 (o'2) as a function of 1 /R for rj = 0.4 and X2 = 0.002 (as given by 
different approximations) and simulation results from Ref. 1221 . Dash-dot-dot line: g1^{cr2)\ dash-dot line: 

922 ("^2); dotted line: gl^i^ {(12); solid line: ^22^ ('''2); dashed line: (;^^((T2); solid circles: simulation data of 
922(1^2) from Ref. (22j. 




of 522 (o'2) in that region seems to be reasonably well captured by both Gf.2+ and 

Ge3-|-- 

Once the question of universality has been examined, in order to complement 
the perspective of the present results, in Figs. ITHlOlwe illustrate some other issues. 
In particular, we want to assess the performance of the theoretical approximations 
obtained with the universality assumption, denoted by 522 ('''2) (where the asterisk 
will be either e2, e3, e2+ or e3+ depending on whether one uses Ge2(^/, z), Gesirj, z), 
Ge2+{j]i ^) oi" G'e3+(?y, -z), respectively, to compute the contact value) and that of the 
AH approximation §22^ {(72) [cf. Eq. ([9])] as compared to simulation data. Hence, in 
Fig. [7] we have plotted 5'22(o'2) as a function of l/R taking X2 = 0.002 and r] = 0.4, 
while in Fig. [8] the large-large contact value has been plotted as a function of X2 



:38 Molecular Physics SYHAH 



10 Santos et al. 




1 

Figure 9. (Colour online) Plot of 522 (o'2) as a function of r) for X2 = 0.002 and ii = 5 (as given by different 
approximations) and simulation results from Ref. [2^ . Dash-dot-dot line: g^{a2)\ dash-dot line: 922 ("^2); 
dotted line: 3%"^^ {a'2)', solid line: g'^^ {02); dashed line: g^^(o'2); solid circles: simulation data of 5122 ("'2 ) 
from Ref. [2^1. 




Q r , I , I , I , I 

0.0 0.1 0.2 0.3 0.4 

Figure 10. (Colour online) Same as in Fig. El but for X2 = 0.002 and R = 10. 

for = 0.4 and R = 5. These plots suggest that the new approximations may 
indeed improve the already reasonable agreement between the AH formula and 
simulation results, but all theories seem to underestimate the height of the peak in 
the 522 (o'2) vs. 1/R plot. The region of interest, namely relatively large R and small 
X2, is further examined in Figs. [9] and [10] which illustrate two additional points. 
On the one hand, one can clearly see that the curves are almost indistinguishable 
for low densities (up to r] ~ 0.2 in these cases) but at higher densities the new 
proposals Ge2+{rj,z) and Ge3+{r],z) indeed improve on the performance of the 
earlier ones Ge2(?7, -z) and Ges^rjjz). On the other hand, they also show that, in 
the high density region, the AH formula performs a little bit better than either 
Ge2-f(??, 2^) or Ge3+{ri,z) as R becomes larger. 
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4. Concluding remarks 

The above results deserve some further consideration. In this paper, following the 
universality approach that some of us have followed for some years, we have in- 
troduced two new proposals [cf. Eqs. ([25]) and (f26l) ] for the contact values of HS 
mixtures which improve on the earlier proposals in the case where one has disparate 
in size mixtures. These proposals fulfill certain exact conditions and are rather flex- 
ible, requiring as their only input the contact value of the one-component system. 



The comparison with recent computer simulation data [22, l29|] indicates that 
the universality assumption may have some limitations but, since there are many 
technical difficulties in simulating mixtures with large R and very small X2, one 
cannot reach definite conclusions on this issue on the basis of these limited data. On 
the other hand, the new proposals following the universality assumptions for the 
large-large contact values seem nevertheless to yield reasonably good agreement 
both with the simulation results and with the recent proposal by Alawneh and 
Henderson [22]. 

It is clear that the AH formulae do not obey the ansatz (jlOp . First, the 
non-universal character of Eq. ([7]) is due to the appearance of the combination 
(a")(cr^)/(cr^)^. In addition, the values given by Eqs. ([8]) and Q depend on 1 — i?^". 
As mentioned above, our analysis of the universality issue in the light of the sim- 
ulation results suggests that the ansatz (jlOp may break down at some point and 
hence whether the AH proposal is universal or not would not be particularly rele- 
vant. However, a more serious problem with the AH formulae is that Eq. ([9|) fails 
to reduce to the exact result to first order in r] given in Eq. ([2]). In contrast, all 
the universal approximations quoted in this paper as well as the VS approximation 
yield such result in the appropriate limit. Whether or not the universality feature 
is confirmed or discarded depends on the availability of further simulation results. 
Our hope is that this paper may encourage the performance of such simulations. 



Acknowledgements 

The work of A. S., S. B. Y. and M. L. H. has been supported by the Ministerio 
de Educacion y Ciencia (Spain) through Grant No. FIS2007-60977 (partially fi- 
nanced by FEDER funds) and by the Junta de Extremadura through Grant No. 
GRU09038. 



References 



[1]J. L. Lebowitz, Phys. Rev. A 133, 895 (1964). 

[2]H. Reiss, H. L. Frisch and J. L. Lebowitz, J. Chem. Phys. 31, 369 (1959); E. Helfand, H. L. Frisch 
and J. L. Lebowitz, J. Chem. Phys. 34, 1037 (1961); J. L. Lebowitz, E. Helfand and E. Praestgaard, 
J. Chem. Phys. 43, 774 (1965). 

[3]M. J. Mandell and H. Reiss, J. Stat. Phys. 13, 113 (1975). 

[4]Y. Rosenfeld, J. Chem. Phys. 89, 4272 (1988). 

[5]M. Keying and D. S. Corti, J. Phys. Chem. B 108, 19756 (2004). 

[6]T. Boublik, J. Chem. Phys. 53, 471 (1970). 

[7]E. W. Grundke and D. Henderson, Mol. Phys. 24, 269 (1972). 

[8]L. L. Lee and D. Levesque, Mol. Phys. 26, 1351 (1973). 

[9]G. A. Mansoori, N. F. Carnahan, K. E. StarUng and J. T. W. Leland, J. Chem. Phys. 54, 1523 (1971). 
[10]D. Henderson, A. Malijevsky, S. Labfk and K. Y. Chan, Mol. Phys. 87, 273 (1996). 
[11] D. H. L. Yau, K.-Y. Chan and D. Henderson, Mol. Phys. 88, 1237 (1996); 91, 1137 (1997). 
[12] D. V. Matyushov and B. M. Ladanyi, J. Chem. Phys. 107, 5815 (1997). 

[13] D. Henderson and K. Y. Chan, J. Chem. Phys. 108, 9946 (1998); Mol. Phys. 94, 253 (1998); 98, 1005 
(2000). 

[14] D. Henderson, D. Boda, K. Y. Chan and D. T. Wasan, Mol. Phys. 95, 131 (1998). 
[15] D. Matyushov, D. Henderson and K.-Y. Chan, Mol. Phys. 96, 1813 (1999). 
[16] D. Cao, K.-Y. Chan, D. Henderson and W. Wang, Mol. Phys. 98, 619 (2000). 



38 Molecular Physics SYHAH 



12 REFERENCES 



[17] C. Barrio and J. R. Solaria, J. Chem. Phys. 113, 10180 (2000). 

[18] D. Viduna and W. R. Smith, Mol. Phys. 100, 2903 (2002); J. Chem. Phys. 117, 1214 (2002). 

[19] D. Henderson, A. Trokhymchuk, L. V. Woodcock and K.-Y. Chan, Mol. Phys. 103, 667 (2005). 

[20]A. Vrabecz and G. Toth, Mol. Phys. 104 1843 (2006). 

[21]A. Takhtoukh and C. Regnant, Fluid Phase Equil. 262, 149 (2007). 

[22] M. Alawneh and D. Henderson, Mol. Phys. 106, 607 (2008). 

[23] Note the typographical error in Eq. (7) of Ref. 22j where a factor of R is missing in the prefactor of 

the exponential correction to the VS formula. 
[24]A. Santos, S. B. Yuste and M. Lopez de Haro, J. Chem. Phys. 117, 5785 (2002). 
[25]A. Santos, S. B. Yuste and M. Lopez de Haro, J. Chem. Phys. 123, 234512 (2005). 
[26] M. Lopez de Haro, S. B. Yuste and A. Santos, Mol. Phys. 104, 3461 (2006). 

[27] M. Lopez de Haro, S. B. Yuste and A. Santos, Alternative Approaches to the Equilibrium Properties of 
Hard-Sphere Liquids, in Theory and Simulation of Hard-Sphere Fluids and Related Systems, Lectures 
Notes in Physics, vol. 753, A. Mulero, ed. (Springer, Berlin, 2008), pp. 183-245. 

[28] D. Henderson, F. F. Abraham and J. A. Barker, Mol. Phys. 31, 1291 (1976). 

[29]A1. Malijevsky, S. B. Yuste, A. Santos and M. Lopez de Haro, Phys. Rev. E 75, 061201 (2007). 



